Longitudinal Flow of Protons from 2-8 AGeV Central Au+Au Collisions 
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Rapidity distributions of protons from central ^^^Au + ^^'^Au collisions measured by the E895 
Collaboration in the energy range from 2 to 8 AGeV at the Brookhaven AGS are presented. Longi- 
tudinal flow parameters derived using a thermal model including collective longitudinal expansion 
are extracted from these distributions. The results show an approximately linear increase in the 
longitudinal flow velocity, {i3^)l, as a function of the logarithm of beam energy. 



PACS numbers: PACS numbers: 25.75.-q, 25.60.Gc, 25.75.Ld 
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Experimental heavy ion programs from the Bevalac, 
SIS, the AGS and the CERN SPS have attempted to 
characterize the distributions of particles emerging from 
high energy collisions in terms of simple thermodynamic 
principles. Testing the underlying assumptions of chem- 
ical and thermal equilibrations at different stages of the 
collision has been attempted by comparing model ex- 
pectations to many different experimental observables. 
Static isotropic thermal emission models applied to ob- 
served particle rapidity distributions from experiments 
at all beam energies consistently fail to describe the ob- 
served shape; such model predictions being too narrow. 
Thermal models which include collective longitudinal ex- 
pansion have been much more successful at reproducing 
the observed rapidity distributions |^ . The additional 
collective motion is attributed to intense pressure gradi- 
ents which develop in the very hot, compressed nuclear 
matter fireballs created in heavy ion collisions. 

At the AGS, rapidity distributions of multiple particle 
species, including pions, kaons, protons and lambda hy- 
perons from central collisions have been simultaneously 
described by a thermal distribution with a common lon- 
gitudinal expansion velocity |^, ^. The agreement be- 
tween the proton and other particle distributions sug- 
gests a high degree of stopping of the incident nuclcons 



at the top AGS energy, which implies even more stopping 
at lower beam energies. However, recent investigations of 
the centrality dependence of the proton rapidity distribu- 
tions from 6-11 AGeV Au-|-Au collisions by E917 suggest 
that the degree of nuclcon stopping may be less than pre- 
viously considered Nevertheless, their flat dN/dy for 
central collisions at all beam energies, fitted by sources 
distributed uniformly in rapidity to y-ycm — 0, could also 
be interpreted in the manner presented herein. 

For the asymptotic case at extremely high beam en- 
ergies, Bjorken proposed |^] that nuclear transparency 
would evacuate the central rapidity region of all of the 
initial nucleons, leaving a hot, high-energy density re- 
gion in which a Quark Gluon Plasma might form. In 
160 AGeV Pb+Pb collisions at the CERN SPS, observed 
net proton ((+)-(-)) rapidity distributions exhibit a 
double-humped character which is consistent with this 
transparency. At SPS energies, the nucleon distributions 
are not describable by a simple thermal model with only 
collective longitudinal flow, but rather need additional 
theoretical consideration of transparency. However, the 
negative hadron (mostly pions) rapidity distributions are 
well-described by this model, which might be interpreted 
as suggesting a significant amount of collective longitu- 
dinal flow for produced particles M . 
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More recently, models which attempt simultaneously 
to characterize many of the freeze-out parameters (trans- 
verse and longitudinal flow, temperature, chemical po- 
tential and source volume) have been attempted on data 
at and above the top energy at the AGS 0. The data 
discussed in this letter are compared only to the locally 
thermal model with longitudinal flow. 

The new data were taken by the E895 Experiment ||] 
at nominally 2, 4, 6, and 8 AGeV (1.85, 3.91, 6.0 and 8.0 
AGeV after correction for energy loss before the target) 
at the Brookhaven AGS using the EOS TPC §. Global 
characterization of the collisions is made possible by the 
nearly 47r center of mass frame coverage provided by this 
detector. Particle identification is achieved by correlat- 
ing the average ionization energy loss in the PIO drift gas 
with the charged particle momenta, reconstructed from 
the curvature of tracks bent by the Multi-Particle Spec- 
trometer (MPS) magnet as they pass through the TPC. 

The 5% most central events were selected using the 
distribution of primary charged particle event multiplic- 
ities. Primary tracks are required to originate within 2.5 
cm of the reconstructed event vertex. This analysis is 
based on approximately 20,000 selected central events at 
each beam energy. 

For particle identification jl^, all tracks are initially 
assigned the proton mass and sorted into rapidity and 
transverse mass bins: 0.1 unit rapidity slices from beam 
to target rapidity, with the mid-rapidity slice correspond- 
ing to \y — Ucml < 0.05, and 25 MeV/c^ wide transverse 
mass bins in the range < mf-mo < 1.0 GeV/c^. The 
mean total momentum and (3j are computed at the cen- 
ter of the (mt-mo,y) bins and are used with a parame- 
terization of the Bethe-Bloch energy loss as a function of 
f3j to fit for the relative particle populations in (dE/dx) 
projected for each (mt-mo,y) bin. The total yield of pro- 
tons in each bin is obtained from the area under the fitted 
proton (dE/dx) distribution. 

Up to a total momentum, p, of 1 GeV/c, the protons 
are isolated in (dE/dx) and easily identified. Above 1 
GeV/c, pions and kaons begin to contaminate the proton 
sample. Between 6 and 8 GeV/c total momentum, the 
proton, deuteron and triton (dE/dx) distributions signif- 
icantly overlap each other and their individual yields can- 
not be resolved. Pion contamination is removed by using 
the observed ratio of oppositely charged pions in the mo- 
mentum range up to 1 GeV/c to extrapolate to higher 
momenta [Q. The relatively cleanly observed negative 
pion yields, combined with the extrapolated pion ratios, 
are used to remove the n'^ contamination from the pro- 
ton distributions above 1 GeV/c. Kaon contamination 
is eliminated by using the measured kaon results from 
0, |l| for the same beam energies. The confusion among 
protons, deuterons and tritons between 6 and 8 GeV/c 
results in a broad hole in the final proton spectra, where 
the proton yields cannot be extracted. Where appropri- 
ate, the reported errors for the raw proton yields include 



an estimate of the additional uncertainty resulting from 
the models used to determine the contamination from 
other particles. 

Simulations of the detector response to protons in all 
regions of phase space were performed using GEANT 
3.21. Small samples (< 4 tracks per event) of proton 
tracks with momentum distributions approximating the 
measured data were embedded into raw events and prop- 
agated through the reconstruction chain. A correction for 
detector efficiency, acceptance and momentum resolution 
was obtained from the ratio of found GEANT tracks to 
the simulated input using the same track quality cuts and 
centrality selection as the data to be corrected. The pro- 
ton corrections were largest at very backward rapidities, 
where detector acceptance causes significant losses, and 
at forward rapidities, where track density is large, and 
at low transverse mass. However, the corrected spectra. 
Fig. |l|, show very good forward/backward rapidity reflec- 
tion symmetry. 

The transverse mass spectra of the protons in each ra- 
pidity slice are shown in Fig. |l|. In order to determine 
yields, these spectra have been fit with the Maxwell- 
Boltzmann distribution. 
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where A{y)mQ represents the mt-mo=0 intercept and 
r(y) is the inverse slope parameter. It agrees fairly well 
at large mt-mo, but the agreement deteriorates for very 
low mj-mg , where transverse flow modifies Eq. (^ . How- 
ever, the resulting systematic overestimate of the total in- 
tegrated yields obtained from the fits are smaller than the 
uncertainties in the measured integrated yields, dN/dy. 

The rapidity density from Eq. (Q) can be obtained by 
integrating over mt with an overall normalization param- 
eter, B, as 



dNth 
dy 



(y) = BT^i-, 
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( — ^ cosh y) 
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(2) 

The width of the distribution, Eq. (^), as a function 
of rapidity for a given emitted particle is a function of 
the temperature and particle mass only. Using the mid- 
rapidity inverse slope parameters (uncertainty of 2 MeV) 
from Eq. ([^) to approximate the temperature (Table |) in 
Eq. (H), this distribution can be compared to the rapidity 
densities extracted from the integrated transverse mass 
spectra. The result, shown as dashed lines in Fig. |[ is 
clearly too narrow to describe our data. 

Schnedermann, et al. modeled this increase in the 
widths as the result of collective longitudinal flow. This 
assumes that the observed distributions arise from the 
superposition of multiple boosted individual isotropic, lo- 
cally thermalized sources at each rapidity, r). This is the 
longitudinally boost-invariant approach |^ but taken to 
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(AGeV) 



{MeV/c^) 



{Pl) mL Ref. 



1.06 (Au) 


80 


0.37 


0.19 


0.194 




1.15 (Au) 


92 


0.41 


0.20 


0.204 1 
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187 


0.570 ± 0.024 


0.28 


0.292 
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211 


0.889 ± 0.024 


0.42 


0.463 
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216 


1.049 ± 0.026 


0.48 


0.547 
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229 


1.104 ± 0.025 


0.50 


0.577 




6 (Au) 


253 


0.990 


0.46 


0.518 , 
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8 (Au) 


267 


1.086 


0.50 


0.577 , 


'41 


10.8 (Au) 


279 


1.166 


0.52 


0.609 , 


41 


11 (Au) 


130 


1.10 


0.50 


0.577 , 
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14.6 (Si+Al) 


120 


1.15 


0.52 


0.609 , 
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158 (Pb) 


160 


1.99 


0.76 


1.169 , 




200 (S) 


220 


1.70 ± 0.30 


0.69 


0.953 , 
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TABLE I: Longitudinal flow parameters extracted from fits of 
Eq. (^ to the proton rapidity densities at 2, 4, 6, and 8 AGeV, 
using the mid-rapidity proton inverse slope parameters, To. 
Note that the parameters extracted at the SPS come from 
produced particles, primarily pions, with comparison among 
other species, such as kaons and strange baryons suggesting 
collectivity. 



a finite boost, rjmax- The distributions are 



dN 
dy 



' dm. . 



(3) 



where rjmax = - flmim from symmetry about the center of 
mass; dNth/dy is Eq. (||) with T from Table | and N is 
the number of protons observed. An average longitudinal 
source velocity can be obtained by averaging over [t^I as 
(/3l) = tanh{r]max/'2)- 

Fits of this form have been applied to a wide array of 
experimental data at the AGS and the SPS §, |, |, 
Although the inverse slope in Table ^ is not the true tem- 
perature, the longitudinal boost, rjmax, has been observed 
to depend only weakly on the temperature input to the 
model, with the exception of nucleons at SPS energies, 
which are affected by transparency [Q. In the present 
data, a 50% decrease in the temperature results in a max- 
imum change of only 15% in rjmax- 

Fits of Eq. (^) to the proton rapidity densities at 2, 4, 
6, and 8 AGeV are also shown on Fig. H The fit parame- 
ter, rjmax, and the corresponding average velocities, (/3l), 
are listed in Table ||. {f3j} l is evaluated by computing 
7 = l/yT^Wf. 

Fig. H shows our results for (Pj) ^ as a function of beam 
energy. Also included are the extracted parameters from 
the compiled data from |, |, g ^ ^ . 

A modification to a previous interpretation of the low- 
est energy FOPI data point is presented here. The pro- 
ton longitudinal flow reported for 1.06 AGeV Au+Au 
in 1^, is from an interpretation of the measured 

data from [|16| which assumed isotropic emission from a 



Siemens and Rasmussen I|T^ radially boosted spherical 
thermal source. Good agreement between their Au-|-Au 
rapidity distributions and predictions from a Siemens- 
Rasmussen style model based on fits to their mid-rapidity 
kinetic energy spectra was reported A temperature 
of 80 MeV and an average (isotropic) flow velocity of (/3r) 
= 0.41 were reported. 

The resulting longitudinal flow velocity of (B) = 0.38 
attributed to the 1.06 AGeV Au-|-Au data by |^ is a fac- 
tor of two larger than the value one extracts using the 
present prescription. The value reported in Table | and 
Fig. H is from our re-analysis of these FOPI 1.06 AGeV 
rapidity densities using their reported temperature of 80 
MeV with Eq. (|). The result is (Pl) = 0.19, which is 
consistent with a value obtained using the same model 
with Au-|-Au collision data measured by the EOS collab- 
oration for E;,eam=l-15 AGeV and reported in |l4[ . 

Based on the reported (/Jl) for Au-|-Au at 1 AGeV and 
11 AGeV, a linear relation between (Pl) and the loga- 
rithm of beam energy was suggested |2| |l^, ^ . Extrapo- 
lated to SPS energies, it did not describe the longitudinal 
flow at the SPS extracted from produced particles such 
as pions. However, with the inclusion of the present data 
from E895 and our re-analysis of the 1 AGeV protons. 
Fig. H shows a linear trend, spanning 1 to 160 AGeV, 
for the heaviest systems, with a slope ~ 2.5 larger than 
was previously concluded. The SPS Pb-|-Pb value (ob- 
tained from produced particles, not protons) no longer 
shows a significant excess above the systematic. Further 
study of rapidity distributions for protons as a function 
of centrality [Q, and produced particles, with system size 
and beam energy in the range above the AGS would help 
to clarify the relation between stopping and longitudinal 
expansion of protons in heavy ion collisions. 
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FIG. 1: Invariant yield per event as a function of mt-mo for protons in central Au+Au collisions at 2, 4, 6, and 8 AGeV. 
Midrapidity is shown unsealed, while the 0.1 unit forward/ backward rapidity slices are scaled down by successive factors of 10. 
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FIG. 2: Proton rapidity distributions in central Au-|-Au col- 
lisions at 2, 4, 6, and 8 AGeV. The dashed curves correspond 
to isotropic emission from a stationary thermal source with 
temperatures given by the mid-rapidity inverse slope parame- 
ters from the transverse mass fits (Eq. (^)), whereas the solid 
curves indicate fits with longitudinal fiow (Eq. (H)). 




FIG. 3: Energy excitation function of longitudinal fiow ve- 
locities ((/37)i), from heavy ion collisions. The open cir- 
cle has been adjusted here as described in the text. For 
heavy systems, a roughly linear dependence with respect to 
\ogifj{Ebeam) ovcr more than two orders of magnitude of beam 
energy is illustrated (solid line) . Note that the parameters ex- 
tracted at the SPS come from produced particles, primarily 
pions, with comparison among other species, such as kaons 
and strange baryons suggesting collectivity. 



